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Abstract 


A time marching Navier-Stokes code called 
PARC3D was used to study the 3-D viscous flow 
associated with an advanced ducted propeller 
(ADP) subsonic inlet at take-off operating 
conditions. At a free stream Mach number of 0.2, 
experimental data for the inlet-with-propeller test 
model indicated that the airflow was attached on the 
cowl windward lip at an angle of attack of 25°, 
became unstable at 29°, and separated at 30°. An 
experimental study with a similar inlet and with no 
propeller (through-flow) indicated that flow 
separation occurred at an angle of attack a few 
degrees below the value observed when the inlet 
was tested with the propeller. This tends to indicate 
that the propeller exerts a favorable effect on the 
inlet performance. During the through-flow 
experiment a stationary blockage device was used to 
successfully simulate the propeller effect on the inlet 
flow field at angles of attack. In the present numer- 
ical study, this flow blockage was modeled via a 
PARC3D computational boundary condition (BC) 
called the screen BC. The principle formulation of 
this BC was based on the one-and-half dimension 


actuator disk theory. This screen BC was applied at 
the inlet propeller face station of the computational 
grid. Numerical results were obtained with and 
without the screen BC. The application of the 
screen BC in this numerical study provided results 
which are similar to the results of past experimental 
efforts in which either the blockage device or the 
propeller was used. 


Nomenclature 

A flow area, m 2 

BC boundary condition 

D diameter, m 

h enthalpy per unit mass, J/kg 

L inlet length, m 

M Mach number 

P total pressure, Pa (N/m 2 ) 

p static pressure, Pa (N/m 2 ) 

prop, propeller 
t static temperature, °K 

W mass flow rate (kg/sec) 
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X distance, m 

a angle of attack, degrees 

S porosity, (A,, r A^J/Aj,, 

A pressure loss coefficient 

circumference position, degrees 
n normal to the surface 

p fluid density, kg/m 3 

Subscripts: 

0 free stream 

blk blockage 

c captured 
eff effective 
fix flux 

loc local value on surface 
pf propeller face 
pk peak 

Introduction 

One of the critical flow characteristics associated 
with a subsonic inlet is the onset of a flow 
separation diming a steep climb at take-off operating 
conditions. Flow separation presents a major 
concern for a subsonic inlet designer. Up to the 
present time, the angle of attack (a) at which the 
inlet flow separation would occur has been generally 
determined by relying on experimental studies. 
Design and development then follow with a 
provision of off-design performance margin of safety 
by increasing the thickness of the inlet cowl lip. 
However, a thick lip would increase the inlet weight 
and at cruise it would develop high drag and result 
in a reduction of the critical flight Mach number for 
the inlet beyond which the divergent drag rise would 
occur. So, it is desirable to keep the inlet cowl lip 
as thin as possible for cruise and as thick as 
necessary to keep the flow attached for take-off. In 
order to achieve an optimized inlet cowl design, 
initially one needs to determine the value of a that 
induces the onset of flow separation associated with 
an existing inlet geometry. Then, the cowl geometry 
can be altered accordingly until the optimized 
thickness is attained. However, problems that arise 
in this design procedure are not simple since 
predicting subsonic inlet flow separation is still a 
challenge in the field of CFD study. 

The CFD community has completed several 
numerical studies of flows associated with through- 


flow and powered nacelles with ducted propellers 
for take-off and cruise operating conditions. Wie et 
al. [1] recently performed an inlet nacelle 
design/analysis study. This study explored the fea- 
sibility of designing a low drag hybrid-laminar-flow 
nacelle by using a combination of an inverse Euler 
calculation and a 2-D flow prediction of laminar- 
turbulent boundary-layer transition. Valarezo [2] 
performed an analysis of a ducted propeller using a 
lower-order velocity panel method to study a 
circulation and pressure loading on propeller blades 
for different gap sizes between blade tip and the 
shroud. Mendenhall and Spangler [3 & 4] 
conducted a numerical study of a ducted fan 
performance by using a potential flow method to 
estimate the aerodynamics at various a’s. Chen et 
al. [5], Uenishi et al. [6], and Hirose et al. [7] 
applied Euler analyses to flow-through and powered 
nacelles including the fan cowl, centerbody, and fan 
exit for axial flow and moderate values of a’s for 
take-off and cruise flight conditions. For take-off 
conditions, Chen analyzed a non-separated flow 
condition associated with a powered nacelle at a 
free stream Mach number (M 0 ) of 0.27 and an 
angle of attack of 25°. Nakahashi [8] developed a 
hybrid method of finite difference (FDM) and finite 
element (FEM) to solve the Reynolds-averaged 
Navier-Stokes (RANS) equations by FDM and 
Euler equations by FEM. This hybrid method was 
applied to predict 3-D non-separated viscous flows 
associated with an inlet nacelle for take-off and 
cruise conditions. The validity of this hybrid 
method was assessed based on the computational 
results only. Srivastava [9] applied a hybrid implicit- 
explicit unsteady scheme to solve the Euler 
equations to allow prediction of the aeroelastic 
characteristics of a ducted propeller at cruise. Hall 
and Delany [10] developed a 3-D Euler-based 
computer code to predict flow characteristics 
associated with a single-stage ducted propeller for 
axial and high a free stream flows. Boldman et al. 
[11] applied a panel method with a built-in com- 
pressibility correction to evaluate several ducted 
propeller inlets for take-off conditions. By 
employing an empirical method based on the inlet 
cowl peak Mach number (M pk ), the results could 
provide an indication of the angle of attack at which 
the inlet flow would separate. The panel method 
did not account for the propeller effect; however, 
without separation the comparison shows that the 
potential flow prediction of inlet static pressures 
agrees favorably with the ducted propeller inlet test 
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data. 


Experimental Background 


Potential flow and Euler analyses are generally fast 
for computer turn-around time solutions. Both are 
capable of providing reasonable flow predictions but 
are limited to certain flow conditions. Potential 
flow calculations are limited to a subsonic non- 
separated flow. Euler calculations can handle both 
subsonic and supersonic flows but are limited to a 
non-separated flow. 

As part of the advanced ducted propeller (ADP) 
research program, NASA Lewis Research Center 
(LeRC) and Pratt and Whitney (P&W) conducted 
experimental studies of the aerodynamics associated 
with several ADP inlet test models with and without 
a propeller. The test data from these studies [11] 
and [12] show that without the propeller the inlet 
flow separated at an a value of a few degrees lower 
than when the inlet was tested with a propeller. 
For a ducted propeller the problems associated with 
the inlet design become more complicated as a 
result of the propeller effect on the inlet internal 
flow performance. The interaction between a 
rotating propeller and a flow stream is very 
complex. As far as a ducted propeller inlet design 
and analysis are concerned, it is important to 
understand the propeller-inlet flow interaction and 
the net results on the overall system aerodynamic 
performance. That understanding could help to 
advance the design and analysis procedures to take 
a full accountability of this propeller-inlet flow 
interaction, thus providing a highly efficient system. 

The present study analyzes an ADP inlet that was 
tested at the NASA LeRC in the 9 x 15-Foot Low 
Speed Wind Tunnel. The study made use of an ex- 
isting general purpose flow solver called the 
PARC3D code which solve the 3-D RANS equa- 
tions. The primary interest of this study was to 
apply the code to predict the ADP inlet 
aerodynamic characteristics at take-off for angles of 
attack that result in three flow conditions which 
include attached flow, unstable flow (just about to 
separate), and separated flow. Grid modeling of 
the ADP inlet did not include the actual propeller 
geometry. The propeller effect in this study was 
accounted for by using a computational BC which 
simulates a blockage effect at the propeller face. 


ADP Simulator - with Propeller 

In 1990, NASA LeRC and P&W jointly conducted 
an experimental research study to develop design 
technology for a low noise/low drag ADP 
propulsion system. An ADP simulator was 
installed in the NASA LeRC’s 9 x 15-Foot Low 
Speed Wind Tunnel as shown in figure 1. The test 
data were taken for a nominal M 0 of 0.2, a’s of 0° to 
35°, and for propeller speeds from 7,500 to 12,000 
RPM. The inlet captured mass flow rates (W c ) 
ranged from 14.06 to 21.12 kg/sec. 

The ADP inlet shown in figure 1 was tested with 
and without the installation of four total pressure 
recovery rakes. The test data shows that the 
removal of these rakes resulted in inlet flow 
separation at an a of about two degrees lower than 
observed when the rakes were installed. For the 
analysis, the computational grid did not model these 
total pressure recovery rakes so it would be 
appropriate to compare the computational results 
with the ADP experimental data that were obtained 
without the installation of these rakes. Figure 2 
shows schematics of the ADP inlet model 
instrumentation without these rakes. Two axial 
rows of static pressure taps are located on the inlet 
cowl windward side (tf = 0°) and on the inlet cowl 
leeward side ( <t> = 180°). A set of pressure taps was 
located circumferentially around the inlet and at a 
location 0.5 cm upstream of the propeller face. 
Instrumentation was not installed on the center 
body because it rotated with the propeller during 
testing. 

Through-Flow and Propeller Simulation 

In 1991 NASA LeRC and P&W conducted an 
experimental study of an inlet without a propeller 
(through-flow nacelle). The effect of the propeller 
on the ADP inlet could be determined by 
comparison of the data from the through-flow 
program with the data from the 1990 test of the 
inlet with a propeller installed. The test techniques 
to simulate this propeller effect were investigated. 
This experimental test program was conducted in 
the United Technologies Research Center’s 10 x 15- 
Foot Subsonic Wind Tunnel. The inlet used in this 
experimental study is called the baseline inlet (Fig. 
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3). The method that was used to simulate the 
propeller effect was to partially block the inlet flow 
by using a stationary blockage device that was 
installed at the propeller face. The device was 
made of 25 tapered rods as shown in figure 4. The 
inlet that was analyzed in the numerical study of 
this paper (Fig. 1) is called the midlength inlet. 
Schematics in figure 5 show a comparison of the 
two inlets. They have the same centerbody and a 
similar cowl geometry with same propeller face 
diameter (D pt ) but the cowl leading edge of the 
baseline inlet is about 5.0 cm longer than that of the 
midlength inlet. The inlet-with-propeller (ADP) 
experimental data reported in [11] shows that 
qualitatively the two inlets have very similar flow 
performance characteristics. 

Figure 6 shows the through-flow static pressure 
distributions on the cowl windward side of this 
baseline inlet at M 0 = 0.2 and W c = 17.41 kg/sec. 
At a = 25°, the inlet flow was attached (Fig. 6(a)) 
and at a = 26° the flow was fully separated (Fig. 
6(b)). The attached flow is identified by a high rate 
of static pressure change or suction around the cowl 
lip producing a low minimum peak static pressure 
at the highlight. The separated flow in this case is 
an inlet lip separation resulting in a relatively flat 
static pressure distribution from the highlight. With 
the inlet lip separated, there is little or no diffusion 
as the flow moves downstream. 

Through-flow data for the baseline inlet shows a 
separation at a = 26°, but test data for the same 
inlet when tested with a propeller indicated that 
separation occurred at a = 30°. This difference 
could be attributed to the effect of the propeller. A 
stationary blockage device (Fig. 5) was installed at 
the propeller face. This device simulated the 
aerodynamic effect of the propeller by partially 
blocking the inlet flow. 

Figure 7 shows a comparison of static pressure 
distributions from the baseline inlet when tested 
with the blockage device and with the propeller for 
M 0 = 0.2 and W c = 17.41 kg/sec. At a = 28° (Fig. 
7(a)) the static pressure distributions are nearly the 
same which indicates that the effect of the blockage 
device was similar to that of the propeller. At a = 
29° (Fig. 7(b)), the static pressure distribution from 
the inlet with the blockage device shows the flow 
was separated while the data for the inlet with the 
propeller shows that the flow was still attached. 


This indicates that the inlet-with-propeller 
performance was at a stall or unstable condition 
from a slight upward shift of the minimum static 
pressure as the angle of attack increased from 28° 
to 29°. At a = 30° both static pressure distributions 
are similar and indicate that the inlet flow was 
separated. This comparison shows that the 
stationary blockage device was capable of simulating 
the effect of the propeller to within one degree of 
inlet separation angle of attack. A comparison with 
the through-flow test results (Fig. 6) shows that 
either the propeller or the stationary blockage 
device was capable of enhancing the inlet 
performance by increasing the flow separation angle 
of attack a few degrees higher than that associated 
with the inlet through-flow. The physics of the inlet 
flow interaction with the propeller or with the 
blockage device is not well understood at this time. 
Physically the stationary rods induced a blockage 
effect on the inlet mainstream flow. Nonetheless, 
the positive contribution of the propeller to the 
angle-of-attack capability could be significant for an 
advanced subsonic inlet design where a compromise 
needs to be made between a thin cowl lip to achieve 
low drag and a thick cowl lip to prevent lip 
separation. This data base suggests that the 
propeller effect can provide a safety margin for a 
ducted propeller inlet at take-off operating 
conditions. 


Computations 

PARC3D Code 

The PARC3D code was selected for this numerical 
study. It is a 3-D, multi-block Reynolds averaged 
Navier-Stokes (RANS) solver. It utilizes the central 
differencing scheme on a generalized curvilinear 
coordinate system. The turbulence model used in 
the code is the Baldwin-Lomax model. An essential 
feature of the PARC3D code is its capability to 
compute flows about complex geometries. 
Computational boundary conditions can be specified 
on any portion of the grid surfaces. The code 
incorporates a semi-automatic time-step control 
function which helps to maintain the stability for the 
flow solution being calculated. Additional details 
about the code can be found in [13], 
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Computational Grid 

A general purpose CFD grid generating tool called 
GRIDGEN [14] was utilized to generate the multi- 
block grid for this study. Figure 8(a) shows a 2-D 
multi-block H-grid, a block of C-grid, and a section 
of the H-grid with embedded C-grid (HC-grid). 
The HC-grid was the grid topology used in this 
numerical study. The reason for employing the 
grid-embedding method was that the H-grid alone 
would induce a computational error in the flow field 
near the inlet hi ghli ght as discussed in [15]. This 
error was eliminated by embedding a C-grid block 
around the inlet cowl. The grid size for each block 
is shown in terms of a 2-D matrix with the first and 
second digits representing the number of grid points 
along the axial and transverse directions, 
respectively. Grid block lines were generated using 
a combination of the Bazier’s curve fit and the 
transfinite interpolation. The internal grid lines and 
clustering at the wall boundaries were generated 
using the algebraic solver followed by an elliptic 
solver to smooth the grid lines as well as to 
rearrange the internal grid spacing. The 3-D grid 
shown in figure 8(b) was made by rotating the 2-D 
grid fifty increments circumferentially through 180°. 
The total number of grid points in the 3-D HC-grid 
is 7.6 x 10 s . 

Computational Boundary Conditions 

Figure 9 shows computational BC’s applied at 
various sections of the grid. The flow field 
throughout the computational domain was viscous. 
The external portion of the flow field terminated at 
the nacelle trailing edge; therefore, the calculation 
did not account for the effect of the plume aft of 
the nacelle. According to an analysis study by 
Uenishi [6], the effect of the plume propagates 
upstream but not far enough to influence the inlet 
internal flow performance. The nacelle surrounding 
flow field was calculated using the computational 
free-stream BC. At the inlet nacelle surface, a no- 
slip BC was specified. The screen BC and the mass 
BC were specified at the inlet propeller face and at 
the nacelle exit, respectively. Mathematical 
formulations of free-stream, no-slip, and mass BC’s 
are reported in [13]. 

Screen BC 


This BC was derived based on the ’one-and-half 


dimensions’ actuator disk theory by Horlock [16]. 
It imposed a flow disturbance that behaves similar 
to a wire-mesh screen. When a flow passes through 
a screen there is some pressure loss. The 
magnitude of this loss depends on the flow velocity 
and the screen porosity. In order to compute the 
flow downstream of the screen, a pressure loss 
coefficient (A) was determined based on an 
experimental study by Pinker and Herbert [17]. The 
fundamental governing equations for the screen BC 
can be written 


Disk 



Pi + /^U] 2 = p 2 + p 2 u 2 2 x-momentum 

p 1 u 1 v 1 = p 2 u 2 v 2 = v, = v 2 y-momentum 

h, + (u 2 2 + v^/2 = h 2 

+ ( u 2 2 + v 2 2 )/2 

h; + u, 2 /2 = h 2 + u 2 2 /2 energy 
Pi/ Ab = P 2 / ^ 2 b state 

A = (pi - p 2 )/(0.5/) 1 u 1 2 ) press, loss coeff. 


The pressure at the exit end of the screen cannot be 
determined analytically and is computed based on 
the value of pressure loss coefficient. The screen 
BC requires a specificaton of the screen porosity 
(5). The value of S is zero for 100 percent flow 
blockage and unity for zero percent flow blockage. 
From [12] based on a one-dimensional flow 
calculation, an amount of flow blockage which was 
effective in simulating the propeller effect induced 
a propeller face Mach number between 0.8 and 1.0. 
Based on this reference, a propeller face Mach 
number of 1.0 was selected to use in a 
determination of S. Using this Mach number and 
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a W c of 21.12 kg/sec would yield a value of 0.72 for 
S. Note that the propeller face Mach number of 1.0 
was chosen simply for simplicity. 

The screen BC is analogous to a passive blockage 
device that does not have any suction capacity to 
generate airflow through the inlet. In order to 
simulate a certain amount of mass flow a mass BC 
is needed downstream of the screen boundary. As 
shown in figure 9, the screen BC was prescribed at 
the propeller face BC6 and a mass BC was at the 
inlet exit BC7. An inlet flow calculation using the 
screen BC at the propeller face means that a mass 
BC at the inlet exit is also required. An inlet flow 
calculation without the screen BC becomes simply 
a through-flow simulation. 

Results and Discussion 

The PARC3D flow calculations were made using a 
computational grid of the ADP midlength inlet (Fig. 
1) for a free stream Mach number (M 0 ) of 0.2, an 
inlet mass flow (WJ of 21.12 kg/sec, and angles of 
attack (a) of 25°, 29°, and 30°. The selection of 
these three values of a was based on the 
experimental test data as shown in [11] which 
indicate that the inlet lip airflow was attached at a 
= 25°, unstable at a = 29°, and separated at a = 
30°. The flow prediction is evaluated by comparing 
the numerical results with the midlength inlet-with- 
propeller (ADP) experimental test data. The effect 
of the propeller was accounted for in this study by 
modeling the propeller blockage effect through the 
use of the computational screen BC specified at the 
propeller face. To generate a desired amount of 
airflow through the inlet, the computational mass 
BC was specified at the inlet exit plane. At a = 25° 
and 29°, the calculations were made with and 
without the screen BC. At a = 30°, the calculation 
was made only with the screen BC. 

Angle of Attack of 25° - No Separation 

Figure 10 shows the inlet flow Mach contours 
obtained from the PARC3D calculations with and 
without the screen BC at the inlet propeller face. 
The Mach contours in this figure for the inlet center 
plane show flow characteristics on the inlet 
windward (lower half) and leeward (upper half) 
sides for M 0 = 0.2, a = 25°, and W c = 21.12 kg/sec. 
Without the screen BC (Fig. 10(a)), the cowl 


stagnation point on the leeward side is at the 
highlight and on the windward side is on the 
external flow surface just downstream of the 
highlight. On the cowl windward side near the 
highlight, the flow sensed a strong suction and 
expanded rapidly to a local peak Mach number 
(MpJ of 1.27 and then diffused downstream. On 
the leeward side near the cowl highlight, the flow 
expanded at a much slower rate than on the 
windward side. The boundary layer thickness on 
the cowl leeward side is thinner than on the cowl 
windward side at locations downstream of the 
highlight. The stagnation point on the center body 
is slightly off the center line; and the boundary layer 
thickness on the upper surface is similar to that on 
the lower surface. With the screen BC, the flow 
field in figure 10(b) shows quantitatively similar 
characteristics to figure 10(a) upstream of the 
propeller face. A comparison of the Mach contours 
near the propeller face on figures 10(a) and 10(b) 
shows the effect of the screen BC on the inlet 
interna] flow field. The use of the screen BC 
resulted in a redistribution of the flow field at the 
propeller face. Therefore, the screen BC apparently 
functioned like a blockage device at the propeller 
face. On the cowl windward side the screen BC 
generated a perturbation in the boundary layer 
which resulted in a reduction of the boundary layer 
thickness immediately upstream and downstream of 
the propeller face. On the cowl leeward side and 
on the center body, the screen BC also generated a 
similar perturbation; but the results are not obvious 
since the boundary layers are very thin. 

Figure 11 shows two sets of circumferential static 
pressure contours obtained from the flow 
calculations with and without the screen BC. Each 
set contains the results for five inlet cross sections 
upstream of the propeller face. A comparison of 
these pressure contours shows the flow blockage 
effect that was induced by the screen BC only 
propagated a very short distance (about 3 cm) 
upstream of the propeller face. Figure 12(a) shows 
circumferential static pressure contours from the 
flow calculation without the screen BC. The levels 
of static pressure lines vary in both radial and 
circumferential directions due to the 3-D flow 
phenomenon as a result of inlet operation at an 
angle of attack. With the screen BC that is shown 
in figure 12(b), the circumferential pressure 
contours display a similar variation. A difference 
between the two figures is near the cowl and the 
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centerbody surfaces on the inlet vertical centerline. 
With the screen BC, there is a large variation in 
static pressure. Without the screen BC, only a 
slight variation in static pressure is apparent. 

The screen BC should affect a ducted inlet flow in 
a m ann er similar to a wire-mesh screen. 
Reference 17 shows the pressure loss that is 
associated with a passage of airflow through a wire- 
mesh screen. The amount of the loss is a function 
of S and the approaching flow velocity. The screen 
BC used in this study was prescribed with a 
constant value of 5. But the calculated flow velocity 
near the propeller face varies as indicated by static 
pressure contours in figure 12. To evaluate the 
characteristics of the screen BC, total pressure 
contours in the streamwise direction were plotted 
for the inlet windward side. These contours plots, 
with and without screen BC, are shown in figure 13. 
A comparison between figure 13(a) with figure 
13(b) shows that the screen BC induced a total 
pressure loss in the flow field downstream of the 
propeller face. This loss varies from about 11 
percent in the mainstream to a small loss in the 
boundary layer flow near the wall. The big loss in 
the mainstream is caused by a high approach flow 
velocity. The small loss near the wall is the result 
of a low boundary layer flow velocity. This evidence 
shows that the screen BC did, in fact, induce a 
similar effect to a wire-mesh screen [17]. 

Figure 14 shows a comparison of inlet cowl static 
pressure distributions between the ADP 
experimental data (with propeller) and 
computational results with and without using the 
screen BC for a = 25°. On the windward side (Fig. 
14(a)), the calculated pressure distribution without 
the screen BC is practically the same as that 
resulting from the use of the screen BC. The two 
predicted static pressure distributions agree 
favorably with the test data. On the leeward side 
(Fig. 14(b)), there is good agreement on the 
external cowl surface, but the predicted distributions 
on the internal cowl surface are slightly lower than 
the test data. The effect of the screen BC on the 
inlet flow, as shown in the Mach and pressure 
contour plots, did not have a significant impact on 
the cowl axial static pressure distribution upstream 
of the propeller face. Figure 15 shows a 
comparison of circumferential static pressure 
distributions on the cowl internal surface at 0.5 cm 
upstream of the propeller face between 


experimental data and computational results. 
Centerbody static pressure distributions with and 
without the screen BC are also shown in figure 15. 
On the cowl side, the calculated distribution without 
the screen BC is higher than when the screen BC 
was used. Without the screen BC the results 
compare more favorably with the test data. The 
difference between the screen and no screen BC 
was due to the flow blockage effect of the screen 
BC. This blockage caused the flow speed to 
increase in order to comply with the required W c 
specified with the mass BC downstream. A local 
increase in flow velocity due to this blockage results 
in a lower static pressure than obtained for the flow 
without any blockage. 

The airflow on the windward lip did not separate at 
a = 25°. Nonetheless, this inlet operating 
condition serves as a test case to study the flow 
characteristics resulting from PARC3D calculations 
of the inlet with and without using the screen BC at 
the propeller face. Evaluation of the computational 
results indicated that the calculations predicted the 
ADP experimental inlet flow reasonably well. Also, 
the results show that the screen BC did simulated a 
flow blockage effect similar to the that associated 
with a wire-mesh screen. However, it has not yet 
been demonstrated how well this computational 
blockage could simulate the propeller effect as 
compared to the effect induced by the stationary 
blockage device used with the baseline inlet (Figs. 
6 & 7). To further assess the effectiveness of the 
screen BC subsequent calculations were carried out 
involving ADP inlet unstable flow (about to 
separate) and flow with separation. 

Angle of Attack of 29° - Unstable Flow 

Figure 16 shows streamwise Mach contours on the 
inlet vertical center plane. The flow characteristics 
for a = 29° are similar in some aspects to those for 
a = 25° that was shown in figure 10. The similarity 
includes locations of stagnation points on cowl 
windward and leeward sides, flow expansion and 
diffusion, and the blockage effect induced by the 
screen BC at the propeller face. The difference 
includes the local peak Mach number and the 
boundary layer thickness particular on the cowl 
windward side. Without the screen BC, the local 
peak Mach number on the cowl windward highlight 
decreased from M plt = 1.27 at a = 25° (Fig. 10(a)) 
to M„ k = 1.25 at a = 29° (Fig. 16(a)). With the 
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screen BC, the peak Mach number increased from 
M pk = 1.29 (Fig. 10(b)) to Mp k = 1.35 (Fig. 16(b)). 
The reduction in the peak Mach number for the 
flow without the screen BC, based on reference 11, 
suggests that the flow at a = 29° may have 
separated. This separation is shown in the velocity 
vector distributions on the inlet cowl windward side 
in figure 17. Both calculations with and without the 
screen BC predicted that inlet-lip flow separation 
which occurred upstream near the highlight and 
continued downstream. With the screen BC, the 
flow reattached about 5 cm upstream of the 
propeller face. Without the screen BC, figure 17(b) 
shows the flow separation extended past the 
propeller face and with the screen BC the flow 
downstream of the reattachment point stayed 
attached. This evidence reveals that the blockage 
effect induced by the screen BC was capable of 
suppressing the separation before it reached the 
propeller face. 

Figure 18 shows circumferential static pressure 
contours for five cross sections at various locations 
upstream of the propeller face. In these plots, the 
effect of blockage induced by the screen BC is 
shown to propagate about 10 cm upstream of the 
propeller face. A comparison with similar static 
pressure contour plots for a = 25° in figure 11 
indicates that the influence of the screen BC for a 
separated flow propagated further upstream than 
that for a non-separated flow. Figure 19 shows 
plots of circumferential static pressure over an inlet 
cross section 0.5 cm upstream of the propeller face 
station. Without the screen BC (Fig. 19(a)), a 
strong circumferential static pressure gradient is 
evident. With the screen BC (Fig. 19(b)), the static 
pressure contours are more uniform in the 
mainstream. In the flow region near the windward 
side, there are localized pockets of pressures that 
depict a complex flow field. This phenomenon 
could have evolved from the flow readjustment from 
separation to reattachment upstream. 

Figure 20 shows a comparison of static pressure 
distributions from computational results with the 
ADP experimental data for M 0 = 0.2, a = 29°, and 
W c = 21.12 kg/sec. On the windward side (Fig. 
20(a)), the pressure distribution for the screen BC 
compares more favorably with the test data than the 
results obtained without the screen BC, particularly 
around the inlet cowl highlight. On the leeward 
side (Fig. 20(B)), similar static pressure 


distributions from the two flow calculations are 
obtained. The distributions compare favorably with 
the test data on the inlet external surface but are 
slightly off on the internal surface. This comparison 
is similar to that for angle of attack of 25°, see 
figure 14(b). Figure 21 shows a comparison of 
circumferential static pressure distributions for a 
location 0.5 cm upstream of the propeller face. On 
the cowl side, the predicted distribution without the 
screen BC is higher than that obtained with the 
screen BC and compares more favorably with the 
test data. On the center body the comparison 
shows a similar trend between the two predicted 
static pressure distributions. The difference is 
similar to the comparison for a = 25° that was 
shown in figure 15. The explanation, as discussed 
earlier for figure 15, was due to the blockage effect 
of the screen BC causing the flow speed to increase 
in order to satisfy the required W c . 

The ADP experimental data indicated that the inlet 
flow separated at a = 30° when the inlet was tested 
with the propeller. The inlet through-flow (without 
any blockage) was expected to separate a few 
degrees earlier. This was substantiated by the 
experimental evidence associated with the baseline 
inlet (Figs. 6 & 7). The PARC3D calculation 
without the screen BC or through-flow simulation 
predicted inlet lip separation which persisted 
downstream past the propeller face. The level of 
the separation cannot be assessed due to the lack of 
the through-flow experimental data for the 
midlength inlet and details of the flow field for this 
inlet. The PARC3D calculation when using the 
screen BC also predicted inlet lip separation, but 
the flow reattached about 5 cm upstream of the 
propeller face where the screen BC was specified. 
The influence from the screen BC blockage effect 
appears to concur with the effect induced by the 
stationary blockage rods used with the baseline inlet 
(Figs. 6 & 7). In order to assess the accuracy of the 
predicted flow field, details data are required. 

Angle of Attack of 30° - With Separation 

Figure 22 shows Mach contours from the PARC3D 
calculation with the screen BC for M 0 = 0.2, a = 
30°, and W c = 21.12 lbs/sec. The flow on the cowl 
windward side stagnated on the outside surface just 
aft of the highlight. From this stagnation point the 
flow accelerated around the lip and reached a local 
peak Mach number of M pk = 1.05. A comparison of 
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this peak Mach number with that for a = 29° (Fig. 
16(b)) shows that it dropped from M pk = 1.35 to 
1.05 as the angle of attack increased from a = 29° 
to 30°, respectively. Downstream of the highlight a 
thick flow boundary layer developed. The screen 
BC induced a blockage effect that tended to force 
the flow to become more uniform at the propeller 
face. It also induced a perturbation in the flow 
boundary layer resulting in a reduction in the 
boundary layer thickness a short distance 
downstream of the propeller face. On the cowl 
leeward side the flow stagnated near the highlight 
and the boundary layer which developed 
downstream on the internal surface is considerably 
thinn er than on the cowl windward side. Figure 
23(a) shows axial velocity vector distributions on the 
inlet cowl windward side. These velocity vectors 
shows that the flow separated just downstream of 
the highlight and continued with no appearance of 
reattachment. Figure 23(b) shows that the flow did 
reattach at a location 1.7 cm upstream of the 
propeller face. This reattachment point is 3.3 cm 
further downstream than for a = 29° as shown in 
figure 17(a). 

Figure 24 shows a comparison of static pressure 
distributions between the ADP experimental results 
and the computational results with the screen BC. 
The experimental static pressure distribution 
indicates that the ADP inlet flow separated. 
Separation is indicated by the relatively flat pressure 
distribution around the highlight. The slight 
increase in pressure downstream of the highlight 
indicates that a small diffusion took place along the 
flow path. With the exception around the highlight 
region, a comparison on the cowl windward side 
shows that the predicted static pressure distribution 
agrees with the experimental data on the external 
and internal surfaces. However, just around the 
highlight the analysis predicted a peak static 
pressure value that is significantly lower than the 
experimental data. The reason for this is not 
readily understood at this point. The comparison 
on the inlet leeward side shows good agreement on 
the external surface and at the highlight but a small 
difference appears on the internal surface. Figure 
25 shows a comparison of the circumferential static 
pressure distributions for a location 0.5 cm 
upstream of the propeller face. The distribution 
from the ADP experimental data is skewed to one 
side indicating that the ADP inlet flow separation 
may be unsteady. The comparison on the inlet 


cowl side shows that the predicted static distribution 
is lower than the experimental data. 

Future Study 

Simulation of Propeller Effect 

The propeller in general is an active device which 
rotates thus generating a centrifugal force on the 
flow and adding energy to the flow. The effects 
induced by the propeller on oncoming flow include 
blockage, swirl, and suction. In this study the 
screen BC only induced a blockage effect on the 
flow field and a perturbation in the flow boundary 
layer. This blockage effect was passive and 
generates neither swirl nor suction on the oncoming 
flow. Instead of adding energy the screen BC 
created a momentum deficit in the flow field. 
Therefore, the blockage effect was the only 
commonality between the propeller and the screen 
BC. The suction in this ADP inlet study was 
generated using a mass BC downstream of the 
screen boundary. So the inlet flow calculation using 
the screen BC still neglected the swirl effect of the 
propeller. This could be a reason that the 
PARC3D calculations with and without the screen 
BC were unable to match the ADP experimental 
pressure distribution on the inlet cowl leeward side. 
A suggestion for future related work would be to 
include the swirl effect into the analysis in addition 
to the blockage effect; e.g., the screen BC. The 
most viable location for this swirl effect is at the 
propeller face. The amount of the swirl to be 
prescribed is left to a future study. 

Prediction of Flow Separation 

Predicting flow separation by CFD is still a 
challenge. The comparison in figure 24(a) shows 
the PARC3D calculation under predicted the inlet 
flow separation. Boundary layer separation has 
been known to evolve from a shock-boundary-layer 
interaction in a supersonic inlet flow, from an over- 
expanding supersonic flow through a nozzle, and 
from a strong adverse pressure gradient in a 
subsonic diffuser flow. With the ADP inlet at an 
angle of attack, the flow around the cowl windward 
highlight involved a rapid expansion from a 
stagnation point followed by a diffusion with a local 
supersonic flow in between. All of these flow 
phenomena occurred over a very short distance. 
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The physics of this type of flow are complex and the 
separation associated with this phenomenon is even 
more complicated. The PARC® code uses the 
Baldwin-Lomax turbulence model which according 
to an analysis done by GJ. Harloff et al. [18] 
predicted a late separation of flow through an S- 
duct. To further evaluate the type of flow 
separation associated with the ADP inlet at take-off, 
a future numerical study is suggested in which a 
different type of turbulence model such as the non- 
equilibrium turbulence closure model (NTM) by 
DA. Johnson and L.S. King [19] and/or the Chen 
and Patel two layer /c-e model [20]. The NTM 
turbulence model requires slightly more 
computational efforts than the Baldwin-Lomax 
algebraic model. The NTM model was developed 
specifically for 2-D flow turbulent boundary layer 
calculations, so it may require a substantial amount 
of work to modify for 3-D flows. 


Conclusion 


A numerical study was implemented to evaluate the 
aerodynamic performance of an ADP subsonic inlet 
at take-off operating conditions with a free stream 
Mach number of 0.2 and angles of attack 25°, 29°, 
and 30°. An existing 3-D RANS flow solver called 
the PARC3D code was utilized in the study. A 
computational mass BC was specified at the inlet 
nacelle exit to induce a required mass flow rate of 
21.12 kg/sec through the inlet. At the propeller 
face a computational screen BC was specified in an 
attempt to simulate a blockage effect induced by the 
propeller. When the screen BC was removed, the 
flow calculation simply represented a through-flow 
inlet (no propeller). The computational results 
were assessed through comparisons with the inlet- 
with-propeller (ADP) inlet experimental data. The 
following findings have been concluded from this 
numerical study; 


indicate that the inlet flow was attached as indicated 
by the ADP experimental test data. The blockage 
effect induced by the screen BC propagated a very 
short distance about 3 cm upstream of the propeller 
face. These results suggest that, for attached flow 
the effect of the blockage downstream may not be 
of any significance for the inlet flow upstream. 

3. At a = 29°, the PARC3D calculations predicted 
flow separation on the inlet lip. The blockage effect 
induced by the screen BC propagated about one- 
half the distance between the propeller face and the 
inlet hightlight, about 10 cm upstream of propeller 
face. With the screen BC, the separated flow was 
forced to reattach at 5.0 cm upstream of the screen. 
Without the screen BC, the separation persisted 
downstream beyond the propeller face location. 
This evidence shows that the blockage effect 
induced by the screen BC was capable of 
suppressing the flow separation after it had 
occurred, thus improving the inlet flow 
performance. These characteristics of the screen 
BC are supported by experimental data for the 
baseline inlet when tested with the blockage rods. 

4. At a = 30°, the numerical results predicted flow 
separation on the inlet lip consistent with the ADP 
experimental test data. Due to the blockage effect 
induced by the screen BC, the separated flow 
reattached at 1.7 cm upstream of the propeller face. 
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Figure 1. Installation of an ADP simulator system in the 
NASA LeRC’s 9 x 15 Foot Low Speed Wind Tunnel. 
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180 ° (leeward) 

Static pressure tap 
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Figure 2. Schematic of the ADP simulator with 
the static pressure instrumentation. 
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Figure 3. Baseline inlet installed in the P&W test facility for 
the experimental study of through-flow and propeller simulation. 
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Tapered blockage 
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Propeller face 
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Figure 4. Blockage rods used at the 
propeller face of the baseline inlet. 
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(b) Midlength; L/D pf = 0.4 


Note: (a) was tested with through-flow and prop, effect 
simulation and (b) was used in this numerical study. 

Figure 5. Comparison of two ADP inlet configurations. 
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(b) After separation, a = 26°. 

Figure 6. Baseline inlet experimental static pressure distributions 
on the cowl windward for through-flow at M 0 = 0.2 and W c = 17.41 kg/sec. 
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2-D H-grid 



H-grid with embedded C-grid 


(a) 2-D grid. 


Figure 8. Computational grids. 
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H-grid with embedded C-grid 

(b) 3-D grid. 


Figure 8. Concluded. 
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Figure 9. PARC BC’s applied in analyses 
using screen BC’s at the propeller face. 
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(b) With screen BC. 
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(a) With no screen BC. 


Effect of screen BC 
at propeller face. 
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Figure 10. Mach contours over the inlet center plane 
for M 0 = 0.2, a = 25°, and W c = 21.12 kg/sec. 
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Figure 11. Inlet circumferential static pressure 
contours at five locations upstream of prop, face 
for M 0 = 0.2, a = 25°, and W c = 21.12 kg/sec. 
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Figure 12. Circumferential static pressure contours 
0.5 centimeter upstream of the propeller face 
for M 0 = 0.2, a = 25°, and W c = 21.12 kg/sec. 
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r - Propeller face 



(a) Without screen BC. 


(b) With screen BC. 


Figure 13. Axial total pressure contours on the inlet 
windward side in region surrounding the propeller face 
for Mq= 0.2, a = 25°, and W c = 21.12 kg/sec. 
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(b) Leeward side, <p = 180° 


Figure 14. Comparison of static pressure distributions 
between analyses and ADP experiment for M 0 = 0.2, 
a = 25°, and W c = 21.12 kg/sec. 
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Static pressure ratio, p loc /P ( 



Figure 15. Comparison of circumferential static pressure 
distributions 0.5 cm upstream of the propeller face for 
M 0 = 0.2, a = 25°, and W c = 21.12 kg/sec. 
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Local supersonic 
M pk = 1.35 

(b) With screen BC. 


Effect of screen BC 
at propeller face. 


Free stream 




Local supersonic 
M pk - 1.25 

(a) With no screen BC. 
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C 0.19 
D 0.22 
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Figure 16. Mach contours over the inlet center plane 
for M 0 = 0.2, a = 29°, and W c = 21.12 kg/sec. 
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(5 cm. upstream of prop, face) 


Without screen BC 



(5 cm. upstream of prop, face) 


With screen BC 

(a) On the cowl windward side. 

Figure 17. Axial velocity vector distributions 
for M 0 = 0.2, a = 29°, and W c = 21.12 kg/sec. 
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(b) Axial velocity vectors around the prop, face 
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Figure 18. Inlet circumferential static pressure 
contours at five locations upstream of prop, face 
for M 0 = 0.2, a = 29°, and W c = 21.12 kg/sec. 
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Figure 19. Circumferential static pressure contours 
0.5 centimeter upstream of the propeller face 
for M 0 = 0.2, a = 29°, and W c = 21.12 kg/sec. 


Page intentionally left blank 




(a) Windward side, 0 = 0° 



(b) Leeward side, 0 = 0° 


Figure 20. Comparison of cowl static pressure distributions 
between analysis and ADP experiment for Mo= 0.2, a = 29°, 
and W c = 21.12 kg/sec. 
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Static pressure ratio, p l0C /P ( 



Figure 21. Comparison of circumferential static pressure 
distributions 0.5 centimeter upstream of propeller face 
for M 0 = 0.2, a = 29°, and W c = 21.12 kg/sec. 
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Figure 22. Mach contours over the inlet center plane 
for M 0 = 0.2, a = 30°, and W c = 21.12 kg/sec. 
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Figure 23. Axial velocity vector distributions 
for M 0 = 0.2, a = 30°, and W c = 21.12 kg/sec. 
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Figure 23. Concluded. 
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(b) Leeward side, <p = 0° 

Figure 24. Comparison of circumferential static pressure 
distributions 0.5 centimeter upstream of propeller face 
for M 0 = 0.2, a = 30°, and W c = 21.12 kg/sec. 
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Static pressure ratio, p loc /P, 



Circumferential position, 0 , deg. 


Figure 25. Comparison of circumferential static pressure 
distributions 0.5 centimeter upstream of propeller face 
for Mq= 0.2, a = 30°, and W c = 21.12 kg/sec. 
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